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ABSTRACT 

Recently correlation analyses between Galactic dust emission templates and a number of 
CMB data sets have led to differing claims on the origin of the Galactic contamination at 
low frequencies, de Oliviera-Costa et al (1999) have presented work based on Tenerife data 
supporting the spinning dust hypothesis. Since the frequency coverage of these data is ideal 
to discriminate spectrally between spinning dust and free-free emission, we used the latest 
version of the Tenerife data, which have lower systematic uncertainty, to study the correla- 
tion in more detail. We found however that the evidence in favor of spinning dust originates 
from a small region at low Galactic latitude where the significance of the correlation itself 
is low and is compromised by systematic effects in the Galactic plane signal. The rest of the 
region was found to be uncorrected. Regions that correlate with higher significance tend to 
have a steeper spectrum, as is expected for free-free emission. Averaging over all correlated 
regions yields dust-correlation coefficients of 180 ±47 and 123 ± 16 /jK /MJy sr^' at 10 
and 15 GHz respectively. These numbers however have large systematic uncertainties that we 
have identified and care should be taken when comparing with results from other experiments. 
We do find evidence for synchrotron emission with spectral index steepening from radio to 
microwave frequencies, but we cannot make conclusive claims about the origin of the dust- 
correlated component based on the spectral index estimates. Data with higher sensitivity are 
required to decide about the significance of the dust-correlation at high Galactic latitudes and 
other Galactic templates, in particular Ha maps, will be necessary for constraining its origin. 

Key words: cosmic microwave background - methods: data analysis - radio continuum: gen- 
eral - radiation mechanisms: thermal and non-thermal 



1 INTRODUCTION 

A small, but significant correlation of existing Cosmic Microwave 
Background (CMB) data with maps of Galactic dust emission has 
been detected. It is important to understand the origin, and hence 
the characteristics of any Galactic emission present in CMB maps 
with structure on angular scales relevant to CMB measurements. 
With such information, we can attempt to remove these contami- 
nating emissions from the data, to be left with a pure CMB map. 

Cross-correlating the COBE DMR maps with DIRBE far- 
infrared maps, Kogut et al (1996a,b) discovered that statistically 
significant correlations did exist at each DMR frequency, but 
that the frequency dependence was inconsistent with vibrational 
dust emission alone and strongly suggestive of additional dust- 
correlated free-free emission (|3yy = —2.15). A number of recent 
microwave observations have also shown that these correlations are 
not strongly dependent on angular scale. Different experiments ( 
DMR: Kogut et al. 1996b; 19GHz: de Oliviera-Costa et al. 1998; 
Saskatoon: de Oliviera-Costa et al. 1997; MAX5: Lim et al. 1996; 
OVRO: Leich et al. 1997), together give, with 95% confidence, 
—3.6 < ^radio < —1-3, consistent with free-free emission over the 
frequency range 15-50 GHz. However if the source of this corre- 



lated emission is indeed free-free, it is expected that a similar corre- 
lation should exist between Ha and dust maps. Several authors have 
found that these data sets are only marginally correlated (McCal- 
lough 1997 and Kogut 1997). Thus, most of the correlated emission 
appears to come from another source. 

Draine and Lazarian (1998a,b) suggest that the correlated 
emission could originate from spinning dust grains. This model 
predicts a microwave emission spectrum that peaks at low mi- 
crowave frequencies, the exact location of the peak depending on 
the size distribution of dust grains, and has a spectral index be- 
tween -3.3 and -4, over the frequency range 15-50 GHz (Draine and 
Lazarian 1998b, Kogut 1999). The spectral indices obtained from 
various experiments do not seem to conform too well to the spin- 
ning dust model alone. However, recently de Oliveira-Costa et al 
(1999) (hereafter DOC99) claim to have found evidence of the spin- 
ning dust origin of DIRBE-correlated emission using the Tenerife 
10 and 15GHz data. They find a rising spectrum from lOGHz to 
15GHz, which is indicative of a spinning dust origin, as opposed 
to a free-free origin, for the dust-correlated Galactic foreground at 
these frequencies. In this paper, we use the most recent Tenerife 
data to show that the correlated emission is not necessarily indica- 
tive of spinning dust. 
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2 THE DATA AND GALACTIC TEMPLATES 

Tenerife is a double differencing experiment that takes data by drift 
scanning in right-ascension at each declination. Data are taken ev- 
ery 1° in RA and at intervals of 2.5° in Dec. We used the Tenerife 
lOGHz and 15GHz data in the region, 32.5° <Dec< 42.5°, and 
0° <RA< 360°, with point sources subtracted and baselines re- 
moved. The Tenerife experiment has beams of FWHM 4.9° and 
5.2° for the 10 and 15GHz experiments respectively, with a differ- 
encing angle of 8.1° between east and west positions on the sky to 
make up the final triple beam (Jones et al (1998)). 

Figure 1 shows the patch of sky covered by the experiment 
in Galactic projection. Figure 2 shows the dust data and the Tener- 
ife data at all the five declinations where data have been taken. We 
have extended the analysis presented in Gutierrez et al (1999) to 
also cover a region outside the high latitude range {b > 40°) that 
they use to constrain CMB fluctuations. Therefore, the data we have 
used are different from that used in DOC99 as they did not extend 
the coverage of the raw data processing to be outside this range. 
This is one of the reasons we would expect differences between 
our results and those of DOC99. A great advantage of using the 
Tenerife data to constrain the spectral index of the dust-correlated 
component lies in the fact that similar angular scales and the same 
patch of the sky can be used. This is not the case when the correla- 
tion from different experiments are compared. Therefore, we have 
not used the two further declinations covered by the 15GHz data set 
(Decs 30° and 45°), but not by the lOGHz data, avoiding any effect 
from features present in these declination strips. However, we have 
redone the calculations using these extra declinations and there was 
no significant difference with the results presented in this paper. 

When looking at Figure 2, a difference between the two sides 
of the RA range, at almost all declinations is apparent. The rise 
(drop because of the double differencing) of the Galactic signal to- 
wards low Galactic latitude seen in the convolved dust maps on 
both sides of the RA range is only followed on one side by the 
Tenerife data. The Tenerife Galactic signal shows a stronger Galac- 
tic centre / anti-centre difference and a variation with declination 
which is different compared to the dust template with a particu- 
larly strong feature at Dec 42.5°. Also note that the noisier lOGHz 
data show a wider spread toward low Galactic latitude, which does 
not appear systematically Galactic, at least when compared to the 
Galactic dust template. This may be due to a systematic error in the 
Tenerife data that is not fully represented by the error bars. 

The Galactic templates used are the destriped DIRBE-l-IRAS 
lOO^um template of dust emission (Schlegel, Finkbeiner and Davis 
1998), the 408MHz survey of synchrotron emission (Haslam et. al. 
1981), and the 1420MHz survey of synchrotron emission (Reich & 
Reich, 1988, hereafter R&R). We used the cleaned and destriped 
versions of the synchrotron emission maps by Lawson, Mayer, Os- 
borne & Parkinson (1987). 



3 CORRELATION ANALYSIS AND RESULTS 
3.1 The method 

Assuming that the microwave data consist of a super-position of 
CMB, noise and Galactic components, we can write 

y = aX+xcMB + n. 

Here is a data vector of A' pixels, X is an A' x M element ma- 
trix containing M foreground templates convolved with the Tener- 
ife beam, a is a vector of size M that represents the levels at which 




Figure 1. The Figure shows the patch of sky covered by the Tenerife exper- 
iment in black, overlayed on the DIRBE template (smoothed to 5°) to show 
the position of the Galactic plane. 

these foreground templates are present in the Tenerife data, n is the 
noise in the data and xcmb is due to the CMB convolved with the 
beam. The noise and CMB are treated as uncorrelated. 
The minimum variance estimate of a is 

1/2 

with errors given by Aa,- = Z,-/ where Z is given as, 

C above is the total covariance matrix (the sum of the noise covari- 
ance matrix and the CMB covariance matrix). The noise covariance 
matrix of the Tenerife data is taken to be diagonal, and the CMB co- 
variance matrix is obtained through Monte-Carlo simulations with 
10000 CMB realizations (this can of course be found analytically as 
well, as in DOC99). In all the above equations X and y are actually 
deviations of the corresponding quantities from the weighted mean 
with weights^ given by C^' (this is different to DOC99 who use 
the actual X and y in the analysis). It should be noted that there may 
be other components of emission present in the data so that the er- 
rors we consider here are only lower limits. Also, the templates that 
we are using may not be perfect and we do not account for possible 
uncorrelated structure in C. If there are other non-Gaussian compo- 
nents present in the data then the minimum variance estimate of a 
will be incorrect. 



3.2 Results 

Listed in Table 1 are the correlation coefficients obtained for the 
Tenerife data sets and the lOO^m DIRBE template, the Haslam, 
and the RScR, with each Galactic template taken individually. Table 
2 shows the correlation coefficients obtained when fits are done for 
two templates jointly. The DIRBE and Haslam correlations corre- 
spond to joint lOO^m-Haslam fits, whereas the R&.R values corre- 
spond to joint lOO/jin — R&R fits. AT = (a±Afl)oca/ is the ampli- 
tude of fluctuations in temperature in the Tenerife data, that results 
from the correlation (Ocal being the rms deviation of the template 
map). The analysis is done for different Galactic cuts. Only posi- 
tive Galactic latitudes are used in these tables as the Tenerife data 
do not extend below a of —32°. 

From Table 1, we see that correlations with each of the three 

* We have tried a variety of different weighting schemes, including uni- 
form weights, to check our method and all results found were within the 
68% confidence limits of those presented in the paper 
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Figure 2. The figure shows the Tenerife and dust data in the RA range that 
lies above the Galactic cut oi h > 20°, for all five declination strips. The 
filled circles and the dashed line represent the Tenerife lOGHz data, the 
empty circles and thin solid line represent the 15GHz data. The data have 
been binned every 5 points in RA. The corresponding lo error bars are 
shown. Units are mK. The thick solid line is the lOO^um data convolved with 
the Tenerife beam. This data has been plotted in units of 0.13 x 1 MJy/sr 



templates are detected with high significance. However these corre- 
lations fall off in significance towards higher Galactic latitudes. For 
comparison, the rms of the lOGHz data (after subtracting noise) for 
b > 20° is 181 ± 9i.iK (the error is due to a 5% calibration uncer- 
tainty). If the templates are not significantly correlated with each 
other, one would conclude that 40% of the signal at 10 GHz is 
Galactic emission from Table 1. It should be noted that the cor- 
relation between the DIRBE and Haslam Galactic templates is sig- 
nificant for all Galactic cuts although if the dust and free-free emis- 
sions (or spinning dust) are 100% correlated then the Galactic emis- 
sion would still only account for 50% of the total emission (assum- 
ing that all of the Galactic foregrounds present at this frequency are 
present in either the Haslam or DIRBE templates). 

We also see that the DIRBE-correlated component has slightly 
higher significance than the Haslam or i?&i?-correlated components 
near the Galactic plane for both the data sets. Away from the Galac- 
tic plane, at 10 GHz, the component correlated to the Haslam map 
becomes dominant (for b > 20°) whereas at 15 GHz the DIRBE 
component remains dominant. However, there is almost no signifi- 
cant correlated emission beyond this Galactic cut. 

3.3 Discussion 

It is seen from Table 1 that individual correlations with dust do 
not show a rising (spinning dust) spectrum with significance be- 
tween 10 and 15 GHz for any Galactic cut. From the joint 100/jm — 
Haslam correlations (Table 2), we see that it appears as if there 
is some evidence for a rising spectrum of dust-correlated emission 
only for the b > 20° cut (as reported by DOC99). Spectral indices 
inferred from the correlation coefficients are listed in Table 3. Er- 
rors are small close to the Galactic plane. However, the values do 
not consistently agree with the spectral indices of any one of the 
three Galactic emission components, -2.8 for synchrotron, -2.1 for 
free-free and positive for spinning dust, for example +1 if the peak 
is at 20 GHz, that could be expected to be present as contaminating 
emissions in CMB data. 

From Table 3 we see that the value of the spectral index in- 
ferred from the DIRBE-correlated component at 10 and 15 GHz is 
negative for most Galactic cuts (the correlations for b > 30° and 
b > 40° are not significant, so it was not possible to infer a spectral 
indexQ). From the positive value of the spectral index for b > 20° 
DOC99 conclude that spinning dust is responsible for producing 
the major part of the DIRBE-correlated emission. However, since 
the spectral indices inferred for all other Galactic cuts are negative, 
a specific model of the spatial distribution of spinning dust would 
be required to explain the data. On the other hand it could also be 
that in the i-range where the Galactic signal in the Tenerife data 
drops the estimation of the spectral index becomes more sensitive 
to systematic effects we have not accounted for. To decide this we 
focus in on the b > 20° range and divide the data further. It should 
be noted that these spectral indices were found using the values 
for a and not AT as the beam sizes between the 10 and 15GHz 
data sets are slightly different giving rise to a lower expected AT 
(although the same relative level of contamination and hence a) at 
15 GHz than at 10 GHz which would systematically bias the values 
obtained. 

When we examine the spectral indices obtained from the com- 
ponent correlated to the Haslam map at the two frequencies, we find 

^ No significant correlations were obtained even for the b > 25° cut, indi- 
cating that only a small region near h 20° Galactic latitude is correlated. 
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that the values from the individual analysis are less negative than 
expected for synchrotron emission (Table 3). Joint correlation of 
the lOO^m dust and Haslam maps with the data result in spectral 
indices that are steeper than expected for synchrotron. This seems 
to imply that most of the emission in the 15GHz data that is cor- 
related with the Haslam map is also correlated with the DIRBE 
maps (this is expected as the region that correlates lies close to 
the Galactic plane), leading to a much lower value of a for the 
synchrotron-correlated component than expected at 15GHz (and 
hence a much steeper spectrum). This may simply be because dust- 
correlated emission is more significant in the 15GHz data. Simi- 
larily since the Haslam -correlated emission is more significant at 
10 GHz, the joint correlation method gives a lower value of a for 
the dust-correlated component and a correspondingly higher value 
for the synchrotron-correlated component. Being able to correct for 
this effect might therefore sort the systematic discrepancy. How- 
ever, since we see this discrepancy with the synchrotron spectral 
index we should be wary of the dust template correlation as well. 

If we now consider the DIRBE template to be an exact pre- 
dictor for either spinning dust or free-free emission (so that all the 
features present in this template are present in the Tenerife data; as 
we are assuming when performing the correlation analysis) then the 
level of 'contamination' a will be independent of the region tested. 
This means that each value of a in Table 1 should be the same 
for the 10 and 15GHz data respectively. We can therefore take the 
variance of a across the different regions analysed as a measure of 
the error in this method, and find that weighted mean and rms at 
10 GHz and 15 GHz are 157 (20) and 122 (10) /;K / (MJy sr"') re- 
spectively. This would correspond to a spectral index of — 0.6i'^Q'5. 

3.4 Further study of the b > 20° cut region 

We now divide the data into halves in RA, which gives two regions, 
one towards the Galactic centre and one towards the Galactic an- 
ticentre, at each declination. Repeating our analysis on these data 
we quantify the visual impression from Figure 2. Although the rise 
of the Galactic plane signal is equally strong in both regions of the 
dust template, we only find a significant correlation with the 10 and 
15GHz data in the Galactic centre region, in spite of the fact that 
the noise errors in both regions are comparable. The results of this 
analysis are listed in table 4. It is seen that of the five declination 
stripes in the centre region only the three at higher declination are 
significantly correlated. The anticentre regions are found to be un- 
correlated or even significantly anticorrelated in two stripes, which 
again demonstrates that the structure is not correlated rather than 
that our sensitivity is reduced when dividing the data. Again if we 
take the weighted mean and rms in the value of a between declina- 
tions, we get 197 (66) and 157 (91) ^/K / (MJy sr"') for the centre 
regions at 10 and 15 GHz respectively, giving a spectral index of 
— 0.6_2 The variation between declinations is high and the spec- 
tral index thus obtained is consistent with both free-free and spin- 
ning dust models. 

We perform yet another split of the data in an attempt to iden- 
tify the region where the correlation mainly comes from. It was 
found that the 20 pixels at lowest Galactic latitude of all decli- 
nations taken together (100 pixels) gave significant correlations, 
while the remaining pixels (750) gave no correlation (see table 5). 
Again the spectral index for dust-correlated emission is small and 
negative, while that for emission correlated to the Haslam map is 
steeper than expected. Since it is clear that the different declinations 
correlate differently, correlating them together is incorrect. Hence, 
shown also in table 5 are the values obtained for the same analy- 



sis on the 20 pixels at lowest Galactic latitude of each declination. 
Even though the signal to noise is much lower here, it is clear that 
the Galactic plane signal correlates with a free-free like spectral in- 
dex while the rest of the pixels shows a significant anti-correlation 
with dust at 10 GHz for some declinations. The Galaxy correlates 
positively while the remaining regions correlate negatively or in- 
significantly. It is clear from analysing such splits that the correla- 
tion coefficients obtained when taking all declinations together as 
in tables 1 & 2 (fc > 20°) and in table 4 & 5, and which seem to 
indicate an almost flat or less negative spectral index between 10 
and 15 GHz, are composed of regions that correlate with a nega- 
tive spectral index and regions that do not correlate significantly at 
all. The anti-correlation, which occurs mostly in the lOGHz data, 
lowers the combined best fit a value at that frequency, and hence 
raises the spectral index of dust-correlated emission. Note here that 
the high values of 'a' obtained need not to be in contradiction to 
the level of free-free allowed by Ha as the detections are all close 
to the Galactic plane, where the level of free-free is expected to 
be high. We still see a large variance between the three declina- 
tions which correlate significantly. If we take the weighted mean 
and rms of the highest three declinations as indicative of a compos- 
ite correlation, we get 313 (172) and 219 (83) /./K / (MJy sr"') 
at 10 GHz and 15 GHz respectively, giving a spectral index of 
—Q.9^\, again consistent with both free-free and spinning dust. 
Note that this value was obtained from the weighted average of 
three pairs of numbers shown in table 5, where each pair indicated 
a negative spectral index. The spectral index obtained simultane- 
ously for the component correlated to the Haslam map is consistent 
within 1 o with the synchrotron spectral index. 

Given that the dust templates correlate with the data only in a 
small part of the Tenerife patch above b = 20° we further test, how 
physically meaningful the detection is. This is important because 
the spectral index of a component in this method is inferred by 
assuming that the difference in the correlation amplitude between 
the frequencies is solely due to the change of emissivity of a sin- 
gle physical component and is not due to a spurious change in the 
strength of the correlation. 



4 SKY ROTATIONS 

Another way to test for systematic errors in the method is to take 
random samples from the Galactic templates themselves as "Monte 
Carlo" simulations. The advantage is that we can test against 
chance correlations with typical structure in the templates without 
having to understand the templates to a degree that would allow 
us to model this structure. Since the rise of the signal towards the 
Galactic plane clearly is a strong feature in the data and high Galac- 
tic latitude correlations are at best weak, we correlate to rotated / 
flipped maps in the northern and southern hemispheres, which all 
have the same orientation relative to the Galactic plane. Figure 3 
shows a and its significance for 10 and 15 GHz, for the b = 20° 
Galactic cut, the Galactic centre region of all declination strips 
taken together, and the spectral indices derived from these values. 
We find that the real patch does not correlate most significantly at 
either frequency. Given the number of patches that correlate more 
significantly we conclude that the real correlation is at most typical 
and only due to the Galactic plane signal. Further we see that the 
15 GHz real correlation although less significant than the neigh- 
bouring points has a higher value of d. It seems that the real dust 
patch, although it does not represent the observed 15 GHz structure 
particularly well, happens to produces a high correlation amplitude. 
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Figure 3. The top four panels show the correlation coefficients and their 
significances when the 10 and 15 GHz data are correlated with different ro- 
tated patches of the dust template. These are joint coirelations of the centre 
regions of all declinations combined, when the Haslam map is held fixed. 
The bottom panel shows the corresponding spectral indices. The first point 
in each panel corresponds to the real patch. 



Indeed the real dust patch shows no strong features (low intensity 
rms), which explains the weak correlation and results in an overes- 
timate of a at 15 GHz. An interesting point is that the mean spec- 
tral index which could be interpreted as the spectral index derived 
from correlations to typical Galactic plane dust signals, correlations 
which are about as significant as the real correlation, has a value of 
— 1.5 ± 1.2, well consistent with free-free emission. The rms vari- 
ation in the value of the correlation coefficient at 10 and 15 GHz 
can be taken as an estimate of the error in a for the real patch. The 
error estimates obtained in this way range from about 20% for a at 
both 10 and 15 GHz for the b>0 Galactic cut, to 64% and 30% 
at 10 and 15 GHz respectively for the Galactic centre regions in 
the b > 20 Galactic cut. As expected, the errors from sky-rotation 
increase with decreasing signal and size of the region and are very 
large for the individual declinations. Note that the lower Galactic 
cuts give the highest correlation with the highest significance for 
the real patch showing that the correlations there are not just due to 
aligned structure of the typical Galactic plane. The error estimates 
presented in this section are similar whether we include or exclude 
the rotated patches near the Galactic centre and anti-centre, where 
significantly different structure could be expected, though in the 
first case errors increase slightly for both 10 and 15 GHz. Hence- 
forth we use results corresponding to the second case. 

Since we have strong indication that the correlations between 
dust and the Tenerife data results only from the Galactic plane sig- 
nal which will be similarly present in many physical components 
of Galactic emission we further investigate this particular feature 
of spatial distribution by modelling it in the next section. 



Dust 15GHz lOGHz Haslam 

Figure 4. The Figure shows the FWHM of the Gaussian model for the 
Galaxy that gave the most significant correlation (or minimum x^) plot- 
ted against the different data sets for the h > 20° Galactic cut and for Dec 
40.0 (triangles) and Dec 37.5 (squares). The model was correlated with the 
Galactic centre regions of the data only. 



5 GALAXY MODELLING 

A problem with the interpretation of the correlation results at dif- 
ferent frequencies purely as a spectral dependence is that it neglects 
the possibility of variations in the shape of the Galaxy due to differ- 
ent components at different frequencies, and the resulting change 
of spatial alignment. For example, this method for obtaining the 
spectral index of correlated emission does not take into account 
the varying extent of the galaxy in the different data sets and tem- 
plate maps, so that the errors that we are quoting in the above ta- 
bles might be too small and could be systematically wrong. If the 
Galaxy has a different spatial extent at 10 GHz, 15 GHz and 3 THz 
(100 /jm), it will result in a higher or lower value for the overall 
correlation coefficient than is actually the case. We investigate this 
effect by modelling the spatial distribution of Galactic emission in 
the next section. 

A significant correlation between the Tenerife data and the 
DIRBE dust template only occurs close to the Galactic plane where 
there is a rise in the overall Galactic signal. This rise may be fre- 
quency dependent if different physical components are present (for 
example free-free and synchrotron) so that care must be taken when 
comparing correlation results between different frequencies. We 
take a simple model of the Galaxy 

/ fc^41n2 
S{b) = So exp =- 

where b is the Galactic latitude, 5 is the amplitude and FWHM is 
the full width half maximum of the Gaussian model for the Galac- 
tic plane. The value of 5o is arbitrary as it will only scale the es- 
timate for a for the various correlation results and not affect the 
relative value between models with different width. Each model is 
then convolved with the appropriate experimental beam, and corre- 
lated with the data, while varying FWHM. 

Figure 4 shows the widths of the Galactic models which cor- 
related best with the data and the Galactic templates for the two 
declination strips in the centre region which do show convincing 
correlation. Two strips showed no significant correlation between 
the data and the templates, and Dec = 42.5° has the strong feature 
at 10 GHz, which is also incompatible with the templates. It can be 
seen that the most significant correlation between the dust/Haslam 
template and the model occurs with Galactic models of different 
sizes (FWHM), which are also different from the sizes obtained for 
the 10 and 15GHz data. However there seems to be a clear trend of 
decreasing width with increasing frequency compatible with more 
synchrotron emission at 10 GHz and more dust-correlated structure 
at 15 GHz. The combined fit to all declination stripes simultane- 
ously also shows this expected trend of increasing Galactic width 
with wavelength. If this systematic change is real we might be able 
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Figure 5. Fitting the structure in the Tenerife data with the Gaussian model and the Galactic templates. Examples are shown for the 10 GHz data at Dec 42.5: 
(a) shows the best fit Gaussian model to the data and residuals (b) shows the Haslam (dashed) and dust (dotted) templates with amplitudes from the joint fit 
result and the thick dot-dashed line is the sum of the two (c) an equally acceptable fit, given the goodness of fit, is achieved when the spectra are normalised 
with the amplitudes of the templates at 15 GHz and are constrained to free-free and synchrotron spectral indices. 




to fit for two distinct components with different width and spec- 
tral index. Although our fits of a single Gaussian component are all 
well acceptable the quality of the data is however not good enough 
to fit for two such components, if they are present independently of 
templates. 

In Figure 4 the points for the dust correlations lie closer to 
the 15 GHz points than to the 10 GHz points, implying a better 
match of the Galactic shape for 15 GHz. When correlating the dust 
template to the data, the significance of any correlation will drop 
and a will change systematically, because the structure at each fre- 
quency is different. When correcting for this mismatch, assuming 
that the correlated component we are looking for has identical ap- 
pearance at any frequency, we find that the correlation amplitudes 
for both 10 and 15 GHz shift upwards. However the 10 GHz am- 
plitude, because of the larger difference in size with dust, gets a 
larger shift, which results in a drop of the spectral index for the 
dust-correlated emission by typically about 0.2. Taking this effect 
of mismatched structure into account for the Haslam correlations 
results in a less steep spectral index, which is more consistent with 
synchrotron emission. 

When we subtract the best fitting Galactic models from the 
data and the templates no correlations remain between the resid- 
uals, neither between the data and the templates nor between the 
templates, showing again that any correlation is only due to the 
Galactic plane signal, which we are able to model. Figure 5 shows 
these best fit models for the 10 and 15GHz data, for Dec 42.5, and 
the residuals. The fit achieved by the single component Gaussian 
model is indeed good and in fact better than the joint correlation fit, 
which we show for comparison. Further we find that since the joint 



fit to the 10 GHz data is rather poor in terms of goodness of fit, 
synchrotron and free-free components with the expected spectral 
indices (-2.8 and -2.1) can be fitted with an almost similar good- 
ness of fit, based on the reduced j^. When the same is done for the 
other declinations the goodness of fit is generally better for both the 
joint fit as well as the model fitting. Further in the other declinations 
also we generally find equally acceptable fits for a dust-correlated 
component which is constrained to give a free-free spectral index. 



6 TEMPLATE ERRORS 

Another yet unaccounted source of error is an intrinsic error in the 
dust template itself, which can be instrument noise or data reduc- 
tion errors due to point source subtraction etc., but also systematic 
variations in the dust emissivity, which might not be followed in 
the same way by the components at 10 and 15 GHz. Note that the 
DIRBE maps at different frequencies are correlated with each other 
to not more than 95%. If we assume that there is a 5% Gaussian er- 
ror in the flux at 100 /jm due to these errors we can model this 
effect to quantify any change that would occur in the correlation. 
To do this we added a 5% error in the flux (on the Tenerife an- 
gular scales) to the maps and calculated the new values of a. We 
performed 300 Monte-Carlo simulations and found that there was 
a systematic drop in a for the 10 and 15 GHz data by factors of 
1.6 and 1.4 respectively, when the centre regions of all declinations 
were taken together for the b > 20° cut (similar values are obtained 
when the entire b > 20° region is taken). Turning this effect around 
for the real data, assuming that there is an intrinsic 5% error present 
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in the template map, we need to increase the values of a given in the 
tables by these factors. The value of a at 10 GHz has to be increased 
relative to 15 GHz, which results in a significantly steeper spectral 
index, again more consistent with free-free emission. However, it 
should be noted that we do not have error maps for the templates at 
present and thus in this current paper we have not attempted a quan- 
titative evaluation of this effect Presumably template errors, given 
that we have an idea of the likely causes, will be greater towards 
the Galactic plane. Since we have found that the only correlations 
we have to report result for pixels close to the Galactic centre, this 
will be an important source of systematic error. 



7 CONCLUSIONS 

Spinning dust emission can be identified and discriminated from 
free-free emission (/indications for a peak in the emission spectrum 
can be found at low microwave frequencies, which are probed by 
the Tenerife experiment. A rising spectrum between 10 and 15 GHz 
can be taken as a prediction of the spinning dust hypothesis, al- 
though the exact location of the emission peak depends on details 
of the model. Tentative evidence for a turnover in the spectrum of 
the Galactic dust-correlated microwave component between 10 and 
15 GHz has been presented by DOC99. We however find that the 
spectral index of dust-correlated emission is negative for all Galac- 
tic cuts except for the b > 20° cut. The Galactic signal, and with it 
the significance of the correlation, decreases with increasing Galac- 
tic latitude, and no correlations are detected in the higher Galactic 
latitude regions of b> 30°. The variance in the value of a between 
regions with different Galactic cuts is rather large. We further find 
that the correlation detected in the b > 20° region comes only from 
a small number of pixels at low Galactic latitude and towards the 
Galactic centre, where signal from the Galactic plane is present. 
This correlation shows a free-free like spectral index, whereas the 
rest of the region was found to be uncorrelated, or even significantly 
anticorrelated. 

Using sky rotations we show that the correlation we see at 
b > 20° is only an alignment of structure due to the rise of the 
Galactic plane signal. Employing a simple model for this struc- 
ture we were able to demonstrate that the spatial distribution of 
Galactic emission is in fact different in the templates and in the 
data, giving rise to a systematic error. We were also able to show 
that this simple model of the Galaxy fits the data generally better 
than the templates. Further, modelling a Galactic free-free compo- 
nent, which correlates with the dust template, generally yields an 
equally acceptable fit to the data. Another significant systematic 
effect arises due to intrinsic errors in the templates, and all these 
effects cause a misleading increase in the inferred spectral index of 
the dust-correlated component between 10 and 15 GHz. 

A comparison of our results with other experiments is pre- 
sented in Figure 7. Here the data points for the Tenerife experi- 
ment, obtained by taking a weighted average of all detections (all 
regions of all Galactic cuts, joint analyses) with errors taken from 
the sky rotations, correspond to values of 180 ±47 and 123 ± 16 
/MJy sr"' at 10 and 15 GHz respectively. Note that the lOGHz 
point on the plot is significantly higher, as compared to that plot- 
ted in DOC99 and the value quoted in our table 2 This is because 
this region consists of parts that are correlated as well as parts that 
do not correlate or even anticorrelate, as in the case of the lOGHz 
data. Here, since we are focussing only on the regions that correlate 
(these regions are the same for both 10 and 15 GHz data and have 
been found to lie close to the Galactic plane), the value is signifi- 




Haslam correlated emission 




10 20 50 100 

Frequency (GHz) 

Figure 7. Summary of our results on the dust-correlated (top panel) 
and synchrotron-correlated (bottom panel) emission compared to results 
from other experiments (COBE DMR-filled ciicles, 19GHz-empty triangle, 
OVRO-filled triangle, Saskatoon-filled square). The straight lines represent 
synchrotron (solid), free-free (dotted) and vibrational dust (dashed) spectral 
indices. The thick solid curve shows a combination of free-free emission 
with a small spinning dust contribution with a peak at 15 GHz in an attempt 
to tit all the data. This tit appears acceptable, as does the single free-free 
component fit, but it should be kept in mind that the data from different 
experiments might not be directly comparable, since they were taken on 
different angular scales and towards different regions of the sky. The Tener- 
ife data points are an average derived from our various results, see text for 
details. 



cantly higher. Note that the high values of a that we get need not 
be in contradiction to the typical level of free-free allowed by Ha 
maps of other regions as our detections are close to the Galactic 
plane, where the level of Galactic emission is expected to be high. 

The spectral index for dust-correlated emission as deduced 
from the 10 and 15 GHz points is less negative (by about la) than 
expected for free-free emission. The spectral index deduced simul- 
taneously from synchrotron-correlated emission is systematically 
steeper than expected. This could be attributed to the effects that 
we have identified and which have all been shown to influence the 
correlation in the same direction. Note also that in this plot only 
the Tenerife and COBE DMR points, at different frequencies, each 
represent data which probe the same angular scales and were taken 
with the same sky coverage. Inferring a spectral index by compar- 
ing different experiments assumes that the Galactic component is 
traced by a given fixed template and does not depend on parameters 
which vary over the sky, but our analysis of the Tenerife data shows 
large variations of Galactic correlation with the sky region. 

With the present analysis of Tenerife data we are not able to 
make a firm claim about the origin of the dust-correlated compo- 
nent, since we do not find convincing support for a spinning dust 
component. This does not have to rule out this hypothesis, since 
environmental conditions or grain sizes, which affect the position 
of the spectral peak, could change systematically with location on 
the sky, particularly in the transition region between low and high 
Galactic latitude. The spatial variation in the correlation amplitude 
and the possibility of the presence of some spinning dust emission 
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along with free-free emission need to be dealt with. The separation 
task is difficult for two reasons. Firstly neither the frequency depen- 
dence nor the spatial distribution of any of the Galactic components 
at low microwave frequencies is particularly well known. And fur- 
ther we would expect these components to be correlated, since we 
find their templates to be correlated, at least at low Galactic lati- 
tudes. The combination of the results from the other experiments, 
shown in Figure 7, might not be strongly constraining, but neverthe- 
less, does not give conclusive evidence for spinning dust emission 
either, without an expectation of the amplitude of free-free emis- 
sion based on dust-correlated Ha emission. 

In order to make more reliable inferences, a pixel by pixel 
separation of components would have to be performed, using the 
Maximum Entropy method for example. In a forthcoming paper 
we shall perform such a separation incorporating information about 
spinning dust. Also, including other data at frequencies lower than 
10 GHz and adding other templates of Galactic emission such as 
the Ha maps from the Wisconsin H- Alpha mapper (Tufte, Reynolds 
and Haffner 1998) would be useful. 
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Table 1. Individual correlations for the Tenerife data, a is the correlation coefficient, and has units fjK{MJy/sr) ' for the 100/jm template, and fjK/K for the 
Haslam and R&.R templates. 



lOGHz 15GHz 



b 


Template 


a it A(a) 


AT^K 


a ± A[d) 


ATfjK 




100^i»i 


214±0 


6435 ±13 


119±0 


3430 ±6 




Has 


591 ±1 


5850±12 


333 ±1 


3189±6 




R&R 


Z0i4U ± o4 


A^Q A _L 1 1 

4oo4 ±11 


141 /O ± J J 


2476 ± 5 


b>5° 




209 ±1 


3127±10 


123 ±0 


1770 ±5 




Has 


545 ±2 


2617 ±9 


323 ±1 


1502 ±4 




R&R 


i /-iZ4 ± oo 


1 nil _L fi 

1 / / j ± y 


lu/yu ±41 


1077 ± 4 


b>lO° 


lOO^im 


268 ±1 


1918±11 


136±1 


967 ±4 




Has 


431 ±3 


1422 ±9 


259 ±1 


839 ±4 




R&R 




IU/9 ± 8 


O C A 1 /IT 

8654 ± 47 


692 ±4 


b>15° 


lOO^HU 


259 ±9 


199 ±7 


132±4 


106 ±3 




Has 


142 ±5 


161±6 


81±3 


90±3 




R&R 


4727 ±148 


207 ±6 


2883 ±84 


125±4 


b>20° 


lOO^im 


131±21 


41±7 


140±13 


43±4 




Has 


71±6 


70±6 


29±4 


28±4 




R&R 


2294 ±256 


54±6 


1514±165 


35±4 


b>30° 


lOOAim 


18±38 


4±8 


32 ±26 


6±5 




Has 


20±8 


20±8 


7±4 


6±4 




R&R 


897 ±390 


17±7 


379 ±255 


7±5 


b>40° 


lOOfjin 


39 ±43 


8±8 


13±30 


2±6 




Has 


18±10 


14±8 


10±7 


8±5 




R&R 


964 ±432 


18±8 


531 ±304 


10 ±5 



Table 2. Joint correlations for the Tenerife data, d is the correlation coefficient, and has units fjK{MJy/sr) ' for the lOOfjm template, and ^jK/K fov the Haslam 
and R&R templates. 



lOGHz 15GHz 



b 


Template 


a±A{a) 


ATfjK 


a±A{a) 


AT^iK 


b>0° 


lOO^jm 


151±1 


4535 ±36 


114±1 


3281 ±19 




Has 


191±3 


1895 ±33 


16±2 


157±18 




R&R 


6829 ±90 


1224 ±16 


1645 ±52 


287 ±9 


b>5'' 


lOOfjm 


154±1 


2311±21 


116±1 


1669 ±10 




Has 


171±4 


821 ±19 


23 ±2 


106±10 




R&R 


5587 ±99 


572 ±10 


1405 ± 54 


140 ±5 


6> 10° 


lOO/jm 


238 ±3 


1704 ±20 


136±1 


961 ±10 




Has 


64±5 


213±17 


0±3 


0±10 




R&R 


1701 ±140 


139±11 


394 ±79 


31±6 


b> 15° 


lOO/jm 


200 ±10 


154±7 


110±5 


88±4 




Has 


94±6 


106 ±7 


26±4 


29±4 




R&R 


3393 ±211 


148 ±9 


1858 ±140 


81±6 


fo>20° 


lOOfjm 


67 ±22 


21±7 


122 ±14 


38±4 




Has 


65±7 


64±6 


12±4 


11±4 




R&R 


1983 ±288 


47±7 


726 ±203 


17±5 


fo>30° 


lOQfjm 


-16±40 


-3 ±8 


18±28 


4±6 




Has 


21±8 


19±7 


6±5 


5±5 




R&R 


1087 ±455 


20±8 


291 ±304 


5±5 


fo>40° 


lOO/jin 


0±49 


0±10 


-21 ±37 


-4±7 




Has 


18±11 


14±9 


13±8 


10±6 




R&R 


1245 ±576 


23 ±10 


842 ±421 


15±8 
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Table 3. Spectral indices obtained from individual and joint congelations. No number is given when the congelation is insignificant. 



Individual analysis joint analysis 

cut dust Haslam dust Haslam 



fo>0° 


_1 4^+0.01 


-i.4i: 


hO.Ol 
-0.01 


-0 69+° '"' 


-6.11; 


H0.34 
-0.37 


fo>5° 


1 T 1+0.01 

-0.01 


-1.29; 


h0.02 
-0.02 


-0 70+'"''* 


-4.95; 


H0.27 
-0.28 


b> 10° 


'■°'-0.03 


-1.25; 


h0.02 
-0.03 


, ,0+0.05 
-l.:SS_(,05 






b> 15° 


-1 66+"'*' 

'■°°-0.16 


-1.38; 


hO.I7 
-0.18 


_ 1 47+0.23 
'•^'-0.24 


-3.17; 


H0.52 
-0.56 


/)>20° 


16+"''^ 


-2.2i; 


hO.54 
-0.56 


1 48+'-3' 
^■^"-i.oo 


-4.17; 


H.OO 
-1.25 


/)>30° 




-2.56; 


h2.35 
-2.95 




-3.09; 


H2.68 
-5.21 


*>40° 




-1.45; 


h3.31 
-4.06 




-o.8o; 


H3.51 
-3.53 



Table 4. Joint conelations obtained for data at different declinations when the data at each dechnation are divided into two halves, one half being towards the 
Galactic centre (marked C) and the other towards the anticentre (marked AC). 



lOGHz 15GHz 



region 


number of pixels 


dust 


Haslam 




dust 


Haslam 


%- 


all Decs, C 


426 


160±33 


89 ±9 


856 


182±22 


27 ±7 


530 


all Decs, AC 


424 


-80±32 


15±10 


481 


43±21 


-50 ±6 


469 


Dec 42.5, C 


88 


225 ± 95 


341 ±27 


259 


261 ±55 


78±18 


122 


Dec 42.5, AC 


87 


-185±73 


29 ±19 


107 


24 ±43 


-12±12 


88 


Dec 40.0, C 


86 


276 ±56 


74 ±33 


94 


225 ±41 


83 ±22 


93 


Dec 40.0, AC 


86 


-148 ±55 


-12±18 


94 


21 ±34 


-1±12 


79 


Dec 37.5, C 


85 


177 ±50 


89 ±22 


100 


173 ±32 


25 ±14 


69 


Dec 37.5, AC 


85 


36 ±52 


37 ±17 


100 


17 ±32 


-5 ±10 


79 


Dec 35.0, C 


84 


75 ±74 


50 ±13 


85 


72 ±45 


12±10 


86 


Dec 35.0, AC 


84 


0±54 


-8±21 


86 


17 ±33 


-13±10 


106 


Dec 32.5, C 


83 


215±94 


-12±20 


57 


-77 ±69 


38±15 


71 


Dec 32.5, AC 


82 


-75 ±92 


7±24 


79 


128±43 


13±11 


89 



Table 5. The following table shows the conelation coefficients a ± A(a) obtained when the 20 pixels with lowest Galactic latitude in the Galactic centre region 
in each declination are correlated with the conesponding structure in the template maps, and those obtained when the remaining pixels are coirelated. All these 
are joint fits. 



lOGHz 15GHz 



region 


number of pixels 


dust 


Haslam 




dust 


Haslam 


x' 


all Decs 


100 


178 ±44 


111 ±14 


479 


167±31 


18±11 


159 


all Decs 


750 


-52 ±27 


22 ±8 


821 


29 ±19 


47±51 


808 


Dec 42.5 


20 


1148 ±374 


279 ±49 


40 


466 ±141 


24 ±29 


24 


Dec 42.5 


155 


- 147 ± 57 


37 ±17 


187 


54 ±35 


-14±11 


174 


Dec 40.0 


20 


264 ±90 


48 ±10 


26 


177 ±59 


35 ±65 


24 


Dec 40.0 


152 


-75 ±47 


0±16 


160 


26 ±30 


10±10 


139 


Dec 37.5 


20 


315 ±144 


163 ±92 


19 


217 ±75 


64 ±54 


10 


Dec 37.5 


150 


10 ±46 


51±15 


180 


6±29 


1±8 


129 


Dec 35.0 


20 


23±121 


53±21 


20 


34 ±77 


3±18 


23 


Dec 35.0 


148 


21 ±48 


14±15 


150 


7±30 


-6±8 


170 


Dec 32.5 


20 


152 ±200 


9 ±58 


18 


32 ±152 


8±42 


10 


Dec 32.5 


145 


34 ±73 


0±17 


125 


79 ±39 


18±9 


155 
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